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W)! Abstract 



We present a rigorous derivation of the flow at arbitrary time in a deterministic 
ceflular automaton model of traffic flow. The derivation employs regularities in preim- 



Q I ages of blocks of zeros, reducing the problem of preimage enumeration to a well known 

lattice path counting problem. Assuming infinite lattice size and random initial con- 
flguration, the flow can be expressed in terms of generalized hypergeometric function. 
We show that the steady state limit agrees with previously published results. 



1. Introduction 

Since the introduction of the Nagel-Schreckenberg (N-S) model in 1992 ||^, cellular automata 
became a well established method of traffic flow modeling. Comparatively low computational 
cost of cellular automata models made it possible to conduct large-scale real-time simulations 
of urban traffic in the city of Duisburg M and Dallas/Forth Worth ITUi. Several simplified 



models have been proposed, including models based on deterministic cellular automata. For 
example, Nagel and Herrmann p considered deterministic version of the N-S model, while 
Fukui and Ishibashi [Q introduced another model (to be referred to as F-I model), which 
can be understood as a generalization of cellular automaton rule 184. Rule 184, one of the 



elementary CA rules investigated by Wolfram |T3[, had been later studied in detail as a 
simple model of surface growth , as well as in the context of density classification problem 
P|. It is one of the only two (symmetric) non-trivial elementary rules conserving the number 



of active sites |I[], and, therefore, can be interpreted as a rule governing dynamics of particles 
(cars). Particles (cars) move to the left if their right neighbor site is empty, and do not move 
if the right neighbor site is occupied, all of them moving simultaneously at each discrete 
time step. Using terminology of lattice stochastic processes, rule 184 can be viewed as a 
discrete-time version of totally asymmetric simple exclusion process. Further generalization 
of the F-I model has been proposed in |^. 

In all traffic models, the main quantity of interest is the average velocity of cars, or the 
average flow, defined as a product of the average velocity and the density of cars. The graph 
of the fiow as a function of density is called a fundamental diagram, and is typically studied 
in the steady state (t -^ oo). For the F-I model, steady-state fundamental diagram can be 
obtained using mean- field argument |^, as well as by statistical mechanical approach ||lT 



or 

by studying the time evolution of inter-car spacing |T2[. In general, little is known about non- 
equilibrium properties of the fiow. In 0], we investigated dynamics of rule 184 and derived 
expression for the fiow at arbitrary time, assuming that the initial configuration (at t = 0) 
was random, using the concept of defects and analyzing the dynamics of their collisions. In 
what follows, we shall generalize results of ||^ for the deterministic F-I traffic fiow model and 
derive the expression for the fiow at arbitrary time. The derivation employs regularities of 
preimages of blocks of zeros, reducing the problem of preimage enumeration to a well known 
combinatorial problem of lattice path counting. Assuming infinite lattice size and random 
initial configuration, the fiow can then be expressed in terms of generalized hypergeometric 
function. We will, unlike in [^, explore regularities of preimages using purely algebraic 
methods, i.e., without resorting to properties of spatiotemporal diagrams and dynamics of 
defects. 

2. Deterministic traffic rules 

Deterministic version of the F-I traffic model is defined on one-dimensional lattice of L sites 
with periodic boundary conditions. Each site is either occupied by a vehicle, or empty. The 
velocity of each vehicle is an integer between and m. If x{i,t) denotes the position of 
the ith car at time t, the position of the next car ahead at time t is x{i + l,t). With this 
notation, the system evolves according to a synchronous rule given by 

x{i,t + 1) = x{i,t) + v{i,t), (1) 

where 

v{i, t) = min {x{i + 1, t) — x{i, t) — 1, m) (2) 

is the velocity of car i at time t. Since g = x{i + l,t) — x{i,t) — 1 is the gap (number of 
empty sites) between cars i and z -|- 1 at time t, one could say that each time step, each car 
advances by g sites to the right ii g < m, and by m sites if ^f > m. When m = 1, this model 
is equivalent to elementary cellular automaton rule 184, for which a number of exact results 
is known P, ||. 

The main quantities of interest in this paper will be the average velocity of cars at time 

t defined as 

1 ^ 
m = j;rJ2^iht), (3) 



and the average flow 0(t) = pv{t), where p = N/L is the density of cars. In what follows, we 
will assume that at t = the cars are randomly distributed on the lattice. When N ^ oo, 
this corresponds to a situation when sites are occupied by a car with probability p, or are 
empty with probability 1 — p. 

In general, if Nk{t) is the number of cars with velocity k, we have 

1 m 

m = j;TT.kN,it). (4) 

fc=i 

When k < m, Nk{t) is just the number of blocks of type 10^1, where O'^ denotes k zeros. 
This means that a probability of an occurrence of the block 10*^1 at time t can be written 
as Pi(lO^l) = Nk/L. Similarly, for k = m, Pi(10™) = Nm(t)/L. As a consequence, equation 
d^) becomes 

-^kmoH)^mmon 
t^i p p 

We will now demonstrate that in the deterministic F-I model with maximum speed m the 

average flow depends only on one block probability. More precisely, we shall prove the 

following: 

Proposition 1. In the deterministic F-I model with the maximum speed m, the average 

flow 4>m(t) is given by 

0„,(t) = l-p-Pt(O"+^). (6) 

To prove this proposition by induction, we first note that for m = 1 equation (j^) gives 
01 (t) = P((10). Using consistency condition for block probabilities Pi(10) + Pt(00) = Pi(0) = 
1 — p, we obtain 0i(t) = 1 — p — P(00), which verifies (^ in the m = 1 case. Now assume 
that (^ is true for some m = n — 1 (where n > 1), and compute (f)n{t)'- 

n-l 

Mt) = nPt(io") + ^jP(ion) = 

= nPt(10") + (n - l)Pt(10"-4) + J2 J^(lOn) 

n-2 

= (n - i)[p,(io"-ii) + Pi(io")] + Pi(io") + Y. jP(ion) 

Using consistency condition P((10"^^l) + Pi(10") = Pi(10"^^) we obtain 

n— 2 

Mt) = Pt(10") + (n- l)P,(10"-i) + ^ jP(lOn) = Pi(10'^) + 0^_i(t) 

i=i 

Taking into account that Pi(10") = Pt(0") — Pt(0"+^) (which, again, is just a consistency 
condition for block probabilities), and using (^ to express (pm-iit), we finally obtain 

0^(t) = l-p-Pi(O™+i). (7) 

This means that validity of (^ for m = n follows from its validity for m = n — 1, concluding 
our proof by induction. 



3. Enumeration of preimages of 0"^+^ 

Proposition 1 reduces the problem of computing (pmit) to the problem of finding the probabil- 
ity of a block of m+ 1 zeros. In order to find this probability, we will now use the fact that the 
deterministic F-I model is equivalent to a cellular automaton defined as follows. Let s{i,t) 
denotes the state of a lattice site i at time t (note that i now labels consecutive lattice sites, 
not consecutive cars), where s{i, t) = 1 for a site occupied by a car and s{i, t) = otherwise. 
We can immediately realize that if a site i is empty at time t, then at time t-|- 1 it can become 
occupied by a car arriving from the left, but not from a site further than i — m. Similarly, 
if a site i is occupied, it will become empty at the next time step only and only if site i + 1 
is empty. Thus, in general, s{i, t + 1) depends on s{i — m,t),s{i — m + l,t), . . . , s{i + 1, t), 
i.e., on the state of m sites to the left, one site to the right, and itself, but not on any other 
site, what can be expressed as 

s{i,t+ 1) = fm[s{i -m,t),s{i - m + 1, t), . . . , s(i + l,t)j, (8) 

where fm is called a local function of the cellular automaton. For the F-I CA, one can write 
explicit formula^ for fm, such as 

fm[s{i — rfiji), s{i — m + l,t), . . . , s{i + 1, t) j = s{i,t) — min{s(i, t), 1 — s{i + 1, t)} 

-|- min I max{s(z — m, t) , s{i — m + 1, t) , . . . , s{i — 1, t)} , 1 — s{i, t) \, (9) 

which, using terminology of cellular automata theory, represents a rule with left radius m 
and right radius 1. In general, after t iteration of this cellular automaton rule, state of a site 
s{i, t) depends on s{i — mt, 0), s{i — mt + 1,0), ... , s{i -|- t, 0), but not on any other sites in 
the initial configuration. Similarly, a block of k sites s{i, t)s{i + l,t) . . . s{i + k) depends only 
on a block s{i — mt, 0),s{i — mt + l,0), . . . ,s{i + k + t,0), as schematically shown in Figure 1. 
We will say that s{i — mt, 0), s{i — mt + 1,0), . . . , s{i + k + t, 0) is an t-step preimage of the 
block s{i,t)s{i + l,t) . . . s{i + k) . Preimages in the F-I cellular automaton have the following 
property: 

Proposition 2. Block 010203 . . . Op is an n-step preimage of a block 0™+^ if and only if 
p = (n + l)(m + 1) and, for every k (1 < k < p) 

k 

E^(«0>o, (10) 

where ^(1) = —m and ^(0) = 1. 

Before we present a proof of this proposition, note that it can be interpreted as follows. 
Let us assume that we have a block of zeros and ones of length p, where p = {n + l){m + l), 
and we want to check if this block is an n-step preimage of a block 0™"^^. We start with a 
"capital" equal to zero. Now we move from the leftmost site to the right, and every time we 
encounter 0, we increase our capital by m. Every time we encounter 1, our capital decreases 

"'^ Since formula (|9|) will not be used in subsequent calculations, we give it withot proof (which is elemen- 
tary). 
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Figure 1: Fragment of a spatiotemporal diagram for the F-I rule with m = 2. States of nine 
sites during three consecutive time steps are are shown, black circles representing occupied 
sites. Block B^ = 101110100 is a 2-step preimage of the block B2 = 100. Outlined sites 
constitute "light cone" of the block B2, meaning that the state of sites belonging to B2 can 
depend only on sites inside the outlined region, but not on sites outside this region. 

by 1. If we can move from Oi to Op and our capital stays always larger than zero, the string 
010203 ... Op is a preimage of 0"^+^. Condition ([T0|) can be also written as 



5Z«i< 



k 



m + 1' 



111 



because ^(x) = 1 — (m + l)x for x G {0, 1}. 
For the purpose of the proof, strings 0102 . 



of length p satisfying ( pTD for a given m 



and for every k < N will be called m- admissible strings. 
Lemma. Let s{l,t)s{2,t) . . . s{p,t) be an m-admissible string. If 

s{i,t + 1) = fm{s{i - m,t),s{i - m + l,t), . . . ,s{i + l,t)j, 



(12) 



and if fm is a local function of the deterministic F-I model with maximum speed m, then 
s{m + 1, t + l)s{m + 2, t + 1) . . . s{p — 1, t + 1) is also an m-admissible string. 

To prove the lemma, it is helpful to employ the fact that the F-I rule conserves the 
number of cars. Let < fc < p and let us consider strings ^i = s(l,t)s(2,t) . . .s{k,t) and 
S2 = s(m + l, t + l)s(2, t+1) . . . s(A;, t + 1). If the string s{l,t)s{2,t) . . . s{k,t) is m-admissible, 
then its first m -\- 1 sites must be zeros. This means that in one time step, no car can enter 
string s(l, t)s{2, t) . . . s{k, t) from the left. On the other hand, in a single time step, only one 
car (or none) can leave the string on the right hand side, i.e.. 



Y^s{i,t)=e+ Y. s{i,t + l) 



(13) 



i=l 



=m+l 



where e G {0, 1}. Three cases can be distinguished: 



(i) All sites s{k — m + 1, t)s{k — 'm + 2,t)... s{k, t) are empty (equal to 0). Then no car 
leaves Si, which means that e = 0, and 

Y.sit,t)= E ^(^,^ + 1)= E ^(^,^ + i)<^-^- (14) 

1=1 i=m+l i=m+l 

The last inequality is a direct consequence of m-admissibility of 5*1. Since the length of the 
string 5*2 is equal to k — m, the above relation (which holds for arbitrary k) proves that 5*2 
is also m-admissible in the case considered. 

(ii) Among sites s{k — m + 1, t)s{k — m + 2,t) . . . s{k, t) there is at least one which is 
occupied (equal to 1), and s{k + l,t) = 1. In this case, since the last site in 5*1 is "blocked" 
by the car at s{k + l,t), again no car can leave string 5*1 in one time step. Therefore, 

k k 

Y.s{i,t)= Yl s{z,t+l). (15) 

i=l i=m+l 



?7i-admissibility of 5*2 implies 



f^ I I k+l k 

>Y^s{t,t) = Y.s{z,t) + l. (16) 

-L ,' 1 ,' 1 



^ + - i=i i=i 



Combining (|T5D with ( 1T6|) we obtain 

sz,t + l <^— , 17 

i=m+l "^ + 1 

which again shows that 5*2 is m-admissible. 

(ill) Among sites s{k — m + 1, t)s{k — m + 2,t) . . . s{k, t) there is at least one which is 
occupied (equal to 1), and s{k + l,t) = 0. In this case, one car will leave right end of the 
string 5*1, therefore 

k k 

Y.s{t,t)= ^ s{t,t + 1)-1. (18) 

i=l i=m+l 

As before, from m-admissibility of Si we have 

iU_|_1 h 7-1 

J:s{^,t) = J:s{^,t)<^, (19) 

^ ^ m + 1 

1=1 1=1 

hence 

E si^,t + l)= J: .(M + l)-l<A±i_i = ^, (20) 

which demonstrates that case (iii) also leads to ?7i-admissibility of 5*2, concluding the proof 
of our lemma. 

Let us now assume that the block Bi = s{l,t)s{2,t) . . . s{p,t) is m admissible (n being 
some fixed integer and p = {m + l){n + 1)) . Applying the lemma to this block we conclude 
that -02 = s{m + 1, t + l)s(2, t + 1) . . . s{p — 1, t + 1) is ?7i-admissible as well. Applying the 
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Figure 2: m-admissible block with rio zeros and rii ones is equivalent to a lattice path from 
the origin to (tt-q, Ui) which does not touch nor cross the line x = my. corresponds to a 
horizontal segment, while 1 to a vertical segment. 



lemma to B2 we obtain m-admissible block B3 = s{2m + 1, t + 2)s(2, t + 2) . . . s{p — 2, t + 2). 
After n applications of the lemma we end up with the conclusion that the string -Bn+i = 
s{nm + 1, n + l)s{nm + 2, n + 1) . . . s{p — n) is ?7i-admissible. Since the length of Bn+i is 
p — n — nm = (n+l)(m+l) — n(m + l) = m + 1, it must, to be m-admissible, be composed of 
all zeros, i.e., B^+i = 0"^"^^. This means that ?7i-admissibility of Bi is a sufficient condition 
for Bi to be an n-step preimage of 0™^^. Reversing steps in the above reasoning, one can 
show that is is also a necessary condition. 



4. Fundamental diagram 

We shall now use proposition 2 to calculate Pt{0"^^^). First of all, we note that Pt{0"^^^) 
is equal to the probability of occurrence of t-step preimage of 0™^^ in the initial (random) 
configuration, that is 

^*(0™+^)=E^o(a), (21) 

where the sum goes over all t-step preimages of 0™"^^. Consider now a string which contains 
no zeros and rii ones. The number of such strings can be immediately obtained if we realize 
that it is equal to the number of lattice paths from the origin to (no, ni) which do not touch 
nor cross the line x = my, as shown in Figure 2. This is a well known combinatorial problem 
0, and the number of aforementioned paths equals 



no — mni /no + Ui 
no + ni I ni ^ 



(22) 



o 
o 
o 




Figure 3: Graph of the probabihty Pt(0™^^) as a function of p for m = 2 and t = 1 (upper 
hne), t = 5 (middle hne), and t = 100 (lower line). 



Probability of occurrence of such a block in a random configurations is, therefore, 

no — mrii friQ + ni^ 



no + ni 



rii 



p"ni-pr, 



(23) 



where p = -P(l). In a t-step preimage of 0™"+^ the minimum number of zeros is 1 + m{t + 1) 
zeros, while the maximum is (n2 + l)(t+ 1) (corresponding to all zeros). Therefore, summing 
over all possible number of zeros i, we obtain 



Pt{0 



m+l^ 



^"+^*+') 2-m[(m + l)(t + l)-^/ (m + l)(t + l) ^ 

~\{m + l){t + l)-i. 



E 

j=l+m(t+l) 



(m + l)(t+l) 
Changing summation index j = i — m{t + 1) we obtain 



Pt{0 



m+l-\ 



(m + l)(t + l) 



E— . 

fr{t + l\ t + l-j 



|p*+W(l_p)-(m)+^-_ 



(24) 



Figure 3. shows a graph of Pt{0"^^^) as a function of p for m = 2 and several values of 
t. We can observe that as t increases, the graph becomes "sharper" at p = 1/3, eventually 
developing singularity (discontinuity in the first derivative) ai p = 1/3. More precisely, one 
can show (see appendix) that 



lim Pi(0 

t— >oo 



m+u ^ f 1 - (m + l)p if p < l/(m + 1), 
I otherwise. 



(25) 



-Poo(0™^^), therefore, can be viewed as the order parameter in a phase transition with critical 
point at p = l/{m+ 1). Using Proposition 1 we can now find the average fiow in the steady 
state 

'^'"^'^^ 1 1-p otherwise, ^^^ 

which agrees with mean-field type calculations reported in as well as with results of 

To verify validity of the result for t < oo, we performed computer simulations using a 
lattice of 10^ sites with periodic boundary conditions. The average fiow has been recorded 
after each iteration up to t = 100 for three values of p: at the critical point p = 1/3 as well 
as below and above the critical point. The resulting plots of the fiow as a function of time 
are presented in Figure 4. Again, the agreement with theoretical curves 

Mt) = 1 - P - E ^ (^T+'i -) 'y'"'-'(^ - pr'"''"' (27) 

is very good. 

Without going into details, we note that the formula (^) can be also expressed in terms 
of generalized hypergeometric function 2F1: 



(1 - p)i+™+™y (1 + m + t + mt)! 

(pm{t) = 1 - p 77- ■ ,w-| , ,.,, : Tvn 2F1 

(1 + m + mt){l + t)l{m + mt)l 



2'-^ .1_1 
2 + m + mt ' p 



{2i 



Since fast numerical algorithms for computing 2F1 exist, this form might be useful for the 
purpose of numerical evaluation of (pmif)- 

5. Conclusion 

We presented derivation of the fiow at arbitrary time in the deterministic F-I cellular au- 
tomaton model of traffic fiow. First, we showed that the fiow can be expressed by the 
probability of occurrence of the block of m + 1 zeros P(0™"'"^). By employing regularities in 
preimages of blocks of zeros, we reduced the problem of preimage enumeration to the lattice 
path counting problem. Finally, we used the number of preimages to find P(0'"''"^), which 
determines the fiow. 

We also found that the fiow in the steady state, obtained by taking t ^ 00 limit, agrees 
with previously reported mean-field type calculations, meaning that in the case of the F-I 
model mean-field approximation gives exact results. This seems to be true not only for the F-I 
model, but also for many other CA rules conserving the number of active sites ( "conservative" 
CA). For example, in |^ we reported that the third order local structure approximation, 
which is a generalization of simple mean-field theory incorporating short-range correlations, 
yields the fundamental diagram for rule 60200 (one of the 4-input "conservative" CA rules) 
in extremely good agreement with computer simulations. Taking this into account, we 
conjecture that the local structure approximation gives exact fundamental diagram for almost 
all "conservative" rules, excluding, perhaps, those rules for which the fundamental diagram 
is not sufficiently "regular" (meaning not piecewise linear). This problem is currently under 
investigation. 
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40 60 100 



Figure 4: Plots of 02(^)/02(oo) as a function of time for p = 0.3, p = 1/3, and p = 0.35 
obtained from computer simulation on a lattice of 10^ sites. Continuous line corresponds to 
the theoretical result obtained using eq. (|28|). 
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Appendix 

In order to find the limit linit^oo PtiO"^^^) we can write eq. p^ ) in the form 

P,(0'"+i) = J2-±-b{t + l-j, (m + l)(t + l),p), (29) 



where 



b{k,n,p) 



.k/ 



p\l-p) 



n—k 



(30) 



is the distribution function of the binomial distribution. Using de Moivre-Laplace limit 
theorem, binomial distribution for large n can be approximated by the normal distribution 



b{k,n,p) 



-{k — npY 
27mp{l-p) ^^^"^ 2np(l - p) ' 



:exp ■ 



(31) 
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To simplify notation, let us define T = t + 1 and M = m + 1. Now, using (|3T| ) to approximate 
b(T — j, MT, p) in (p9D , and approximating sum by an integral, we obtain 

P^ 0^+M = / = exp ^ ^dx. 32 

^ ^ A T ^2nMTpil - p) 2MTp(l-p) ^ ^ 

Integration yields 



Mp(l-p)r /-(l-T + MpT)2\ ( -MpT\ 



PtiO-^^l = \ ^rb^ exp \,,^ , "/ - exp --^ + 



2nT 



2MTp(l-p) j """P 1,2(1 -p)j 

^(1 - MP) (erf f^ ^^^ ) - erf f J_£±££L 

2 [ VJ2Mp(l-p)Ty VJ2Mp(l-p)T, 



where erf (a;) denotes the error function 

2 f^ 



erf(x) = ^ / e-* dt. (33) 

vr JO 



The first term in the above equation (involving two exponentials) tends to with T —>■ oo. 
Moreover, since lim^.^ooerf(x) = 1, we obtain 

lim Pt(0™+M = i(l - Mp) J 1 - lim erf ( J_^£±^££ 
*-oo *^ ^ 2^ ^M T^- \./2Mp(l-p)T, 



Now, noting that 



„,„erfl4^I±£C =|\ 'ff'-Sl- (34) 

T^oo \ /2Mp(l - p)TJ I -1' otherwise. 



and returning to the original notation, we recover eq. (E 



lim P fO'"+M = I ^ ~ ^"^ ^ ^^^ ifp<l/(m + l), 

t^oo *^ '^ I otherwise. ^ ^ 



References 



[1] Nino Boccara and Henryk Fuks. Cellular automaton rules conserving the number of 
active sites. J. Phys. A: Math. Gen., 31:6007-6018, 1998, |adap-org/9712003| . 



[2] J. Esser and M. Schreckenberg. Microscopic simulation of urban traffic based on cellular 
automata. Int. J. Mod. Phys. C, 8:1025-1036, 1997. 

[3] Henryk Fuks. Solution of the density classification problem with two cellular automata 
rules. Phys. Rev. E, 55:2081R-2084R, 1997, |comp-gas/ 97051101. 



11 



[4] Henryk Fuks and Nino Boccara. Generalized deterministic traffic rules. Int. J. Mod. 
Phys. C, 9:1-12, 1998, |adap-orK/9705003 . 



[5] M. Fukui and Y. Ishibashi. Traffic flow in ID cellular automaton model including cars 
moving with high speed. J. Phys. Soc. Japan, 65:1868-1870, 1996. 

[6] J. Krug and H. Spohn. Universality classes for deterministic surface growth. Phys. Rev. 
A, 38:4271-4283, 1988. 

[7] Sri Gopal Mohanty. Lattice Path Counting and Applications. Academic Press, New 
York, 1979. 

[8] K. Nagel and M. Schreckenberg. A cellular automaton model for freeway traffic. J. 
Physique /, 2:2221-2229, 1992. 

[9] Kai Nagel and Hans J. Herrmann. Deterministic models for traffic jams. Physica A, 
199:254-269, 1993. 

[10] P. M. Simon and K. Nagel. Simplified cellular automaton model for city traffic. Phys. 
Rev. E, 58:1286-1295, 1998, |cond-mat/ 9801022 . 



[11] Bing-Hong Wang, Yvonne- Roamy Kwong, and Pak-Ming Hui. Statistical mechanical 
approach to Fukui-Ishibashi traffic flow models. Phys. Rev. E, 57:2568-2573, 1998. 

[12] Bing-Hong Wang, Lei Wang, P. M. Hui, and Bambi Hu. Analytical results for steady 
state of traffic flow models with stochastic delay. Phys. Rev. E, 58:2876-2882, 1998, 
|cond-mat/ 9804269 . 



[13] S. Wolfram. Cellular Automata and Complexity: Collected Papers. Addison- Wesley, 
Reading, Mass., 1994. 



12 



